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The metallicity and dust-to-metals ratio of the Galaxy are fundamental parameters in understanding the interstellar medium (ISM). 
Currently, there is still some uncertainty surrounding these parameters. In this paper, the dust-to-metals ratio in the Galaxy is deter- 
mined using the photoelectric absorption of the soft X-ray afterglows of a large sample of several hundred gamma- ray bursts (GRBs) 
to determine the metal column density in combination with Galactic dust maps to determine the line-of-sight dust extinction through 
the Galaxy in the direction of the GRB. GRB afterglows often have large extragalactic soft X-ray absorptions and therefore the GRB 
sample's upper-bound will define the Galactic dust-to-metals relation. Using a two-dimensional two-sample Kolmogorov-Smimoff 
test, we determine this upper-bound and so derive the dust-to-metals ratio of the Galaxy. We find Ay = 2.2*^-^ x 10^' cm"- A'h as- 
suming solar, Anders & Grevesse (1989), metallicity. This result is consistent with previous findings using bright X-ray sources in 
the Galaxy. Using the same technique but substituting the H i maps from the Leiden-Argentine-Bonn survey for the dust maps, allows 
us to place a limit on the metallicity in the Galaxy. We find a metallicity consistent with the Anders & Grevesse (1989) solar values 
often used in X-ray fitting. Based on this and previous studies, we suggest that the metallicity of a typical ISM sightline through the 
Galaxy is ~ 0.25 dex higher than the current best estimate of the solar metallicity. We further show that the dust-to-gas ratio seems 
to be correlated with the total gas column density, and that this may be due to the metallicity gradient observed toward the Galactic 



centre. Based on the non-constant nature of the dust-to-gas ratio, we propose that the dust column density, at A'h ■ 
represents a better proxy for the soft X-ray absorption column density than Hi maps. 
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1. Introduction 

The relationship between gas, metals and dust that defines the 
interstellar medium (ISM), plays a central role in the properties 
of star-formation, and in the appearance, evolution and ultimate 
fate of galaxies. However the basic quantitative relationship be- 
tween gas, metals and dust is still not well-defined even for our 
own galaxy. Many studies have been made over the years ex- 
amining the dust-to-gas ratio in the Galaxy and more recently at 
cosmological distances. Two methods dominate these analyses: 
a) comparing hydrogen absorption to dust extinction or redden- 
ing, using the Lye line in the UV (and sometimes H2 UV lines) 
to get the total gas column and pairs of stars to obtain the total 
extinction, and b) soft X-ray photoelectric absorption that mea- 
sures the total metal column density in the foreground of bright 
X-ray sources, and converting this to a gas column assuming 
a metallicity, and comparing this to a dust extinction obtained 
from methods like the Balmer decrement, the deviation in the 
optical/infrared from a blackbody spectrum or even measuring 
the dust column using the halo made by small angle scattering 
of X-rays off the dust. The first method provides a genuine gas- 
to-dust ratio, in the sense that it measures the bulk of the gas 
directly. However, it suffers from the deficiency that it is insensi- 
tive to ionised gas and unless the molecular lines are measured, 
also to H2, where the fraction of hydrogen in H2 may be ~ 0.5 
for Ay > 0.5 ( |Rachford et al.|2009[ ). Furthermore it only works 
along relatively low extinction lines of sight (Ay < 2 - 3) since 
it requires spectroscopy in the UV where the extinction is far 
higher and stars are often too faint to observe in Lya or H2. 



The X-ray absorption method works to very high extinctions 
(Ay ~ 30 or more) and measures essentially all metals whether 
they are ionised or even in the solid phase, providing a census 
of the total column density in metals. However, it is a measure- 
ment of the metal column, not the gas column since the X-ray 
absorption is almost insensitive to hydrogen absorption and de- 
pends only weakly on helium. Therefore the X-ray absorption 
requires a metallicity conversion to move from a dust-to-metals 
to a dust-to-gas ratio. 

It is perhaps worth noting that most studies to date in the 
Galaxy have either focussed on relatively nearby sets of ob- 
jects or provided a small number (between 3 and about 20 for 
the X-ray studies) of lines of sight to locations within a few 
kpc of the Sun. Most studies in the Galaxy have therefore not 
probed the ISM of the Milky Way in either a complete or un- 
biased way. These studies have consistently found a dust-to-gas 
ratio at a level of Nh/ Ay ~ 2 x 10^^ cm~^mag~^ ( BohUn et al. 



1978; Whittet 1981; Diplas & Savage 1994; Ryter 1996; Rein_a 
& Tarenghi 1973; Gorenstein 197 5; Ryter et al. 1975; Predehl & 
Schmitt 1995; Vuong et al. 2003 ; jGiiver & Ozel 2009, Rachford 
et al.|2009_) . The statistical errors quoted for some studies have 
been as small as a few lO'^cm^^mag"', however the variation 
from study to study is closer to a few parts in 10^** cm"^ mag"' 
This discrepancy may be related to an underestimate of the un- 
certainties or to an intrinsic scatter in the relation. 

Gamma-ray bursts (GRBs), while found at cosmological dis- 
tances, are extremely bright. In this paper we use the large ho- 
mogeneous sample of Swift GRB X-ray afterglows to determine 
upper limits to the metal column densities to several hundred 
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lines of sight through the Galaxy. We compare these metal col- 
umn densities to all-sky hydrogen and dust surveys to obtain a 
new dust-to-metals ratio and metallicity value for the Galaxy. 
The GRB afterglows are subject to absorption by their hosts, 
with a minor contribution from intervening objects, so we use 
a 2D 2-sided Kolmogorov-Smirnov (KS) test to overcome this 
limitation to define the Galactic lower envelope to their absorb- 
ing column densities. Since it is X-ray selected, the sample is not 
afflicted by observational bias as UV studies are. However, its 
greatest benefit over previous samples is that this sample passes 
lines of sight at random through the Galaxy, and passes through 
the entire Galaxy in every direction, providing a set of sightlines 
less affected by the relatively local nature of some previous stud- 
ies. 

In the next section I describe the sample, data reduction and 
analysis techniques. In section [3] I present the results of the 
study. Section |4] contains an analysis of the relevance of the re- 
sults and a comparison to previous efforts inside and outside our 
Galaxy. In section [5] I offer my conclusions. All errors quoted 
are statistical uncertainties at the 68% confidence level for one 
parameter of interest unless stated otherwise. 

2. Sample selection, data reduction and analysis 

The aim of the work is to obtain equivalent hydrogen column 
densities (in essence the metal column density, A^Hx) for a large 
number of sightlines through the Galaxy from the soft X-ray 
photoelectric absorption of GRB afterglows. These column den- 
sities will be a combination of the whole column density along 
the line of sight to the GRB, dominated by the Galactic column 
and the absorption from the GRB's host galaxy. This would yield 
a distribution of column densities with a lower limit at the val- 
ues of the Galactic column density. By obtaining the Galactic 
dust and H i column densities (E(B - V)) and A^hi respectively) 
along the line-of-sight to the GRB from surveys, specifically 
SchlegeletaT] ( [T998l l for Galactic dust and |Kalberla"etaL] ( |2005] l 
for H I, the dust-to-metals ratio and metallicity can be obtained 
by suitably fitting the two-dimensional distributions of Nux- Ay 
and Nhx-Nhi with a cut-off' at the value of the Galactic relation. 
An Rv - 3.1 was assumed to convert E{B - V) to Ay, consis- 
tent with the mean and median Ry values found by Fitzpatrick 



found in the conversion, i.e. in Ry, is 0.27 ( [Fitzpatrick & Massa] 

|2007l ). 

The 21cm maps ([Kalberla et al. 2005 [l are the combined 



Leiden/Dwingeloo (Hartmann & Burton 1997[l a nd Instituto 
Argentino de Radioastronoma ( |Bajaja et al.||2005[ ) surveys de- 
tecting H I over a velocity range of -450 km s~^ to -1-400 km s"' 
at 1 .3 km s"' resolution with very high equivalent main beam ef- 
ficiency (> 0.99). 

To obtain an equivalent hydrogen column density, Nux, from 
low-resolution X-ray spectroscopy, the deviation from a power- 
law at low energies is measured. This deviation is due to pho- 
toelectric absorption by metals (primarily O, C, Si, Fe and He 
depending on the energy), and is fit with absorption models of 
the gas. I assume a neutral gas, although the total column density 
is not strongly affected by the first few ionisations of the metals 
since the absorption is due to inner shell electrons. Measured 
cross-sections are used for the elements, and while improve- 
ments to the cross-sections over the years have improved fits to 
absorption data ([Wilms etal.|2000| WAMOO), the effects on low- 



resolution spectroscopy are not large, a few percent at most on 
the total column density (WAMOO). The most important effect 
is the assumed abundance of the elements relative to hydrogen 
and this is discussed in section |3]below. The relative abundances 
of the various metals also has some effect on the total derived 
column density, however, at these resolutions, there is no way to 
distinguish which elements dominate the absorption, and since 
the majority of these elements contribute to the dust, it is not a 
crucial point. 

Recent work suggests that the hot intergalactic medium 
(IGM) might increase the X-ray absorption smoothly with red- 
shift at a low level ( [Behar et al.|2011) . However 1) any such 
effect is small, 2) does not appear to be consistent with IGM ab- 
sorption since no excess absorption is detected in about a third 
of quasars examined as a comparison sample, and 3) may be re- 
lated to inadequately modelled Galactic absorption. Apart from 
smooth IGM absorption, intervening galaxies are unlikely to 
contribute much to the X-ray absorption in the general case as 
the observed absorption is related to the total metal column den- 
sity and drops approximately as (1 -i-z)^-'' due to bandpass effects. 
Every line of sight would have to have the equivalent of a 0. IZq 
metallicity absorber at redshift z = 1 with a neutral gas column 
&Massa [ (2007 [) for 328 sources to distances up to ~ 5 kpc dis- density of ~ 5 x 10^' cm"^ in order to begin to impact the result 

here. Such large column density intervening systems do not exist 
along most sightlines, as argued in Watson et al. (200 7 ). 

All GRBs observed with Swift's X-ray telescope (XRT) 
were included in the sample up to the end of November 2010 
(GRB 101030A), resulting in 638 GRBs. While X-rays are de- 
tected for almost all GRBs observed by Swift-XRT, a significant 
fraction do not have sufficient signal-to-noise to provide a spec- 
trum good enough to determine the X-ray absorption. 

The method adopted here has been to take all pre-reduced 
spectra from the Swift/XRT GRB spectrum repository (Evans] 



tributed across the Galactic plane. 

The dust column values were obtained from Schlegel et al. 
( 1998 1 and reduced by 14% following the analysis of [Schlafl^ 
[et al.[ ( |20I0| l, with uncertainties based on the standard devia- 
tion of nearby points. The atomic hydrogen column densities 
were obtained from the Leiden- Argentine-Bonn survey reported 
in [Kalberla et al.[ (2005) and accessed through the nH tool in 
Ftools. Uncertainties in the Hi column density were set at 10% 
( [Wakker et al.T20T0l i. 

To derive reddenings, [Schlegel et al.| ( [I998| l use the lOO/zm 
sky emission maps from COBE/DIRBE and IRAS/ISSA with 
temperature corrections. The conversion from far-infrared (FIR) 
flux to reddening is made using the excess colours of a sample 
of elliptical galaxies. The zodiacal light contribution is removed 
using the DIRBE 25 //m maps. The values of E{B - V) reported 
should be fairly accurate, since the conversion from FIR emis- 
sion is measured on reddenings, with systematic uncertainties at 



et al. 2009[) for both the windowed timing (WT) and photon 



a level below our ultimate statistical uncertainty ( Schlafly et al 



[2010 ). The accuracy of using a conversion from reddenings to 
absolute extinction in the Galaxy is a matter of significant debate 
(e.g. Fitzpatrick & Massa|2007 1, and will evidently be different 
along different lines of sight. However, the standard deviation 



counting (PC) modes using the appropriate corresponding cal 
ibration files from that archive. These data were then fit with 
a power-law with photo-electric absorption (phabs(pow) in 
Xspec). 

While the PC and WT mode spectra for a given GRB are 
not independent objects, they are independent measurements of 
the same sightline. It is well known that the absorption in the X- 
ray afterglows of some GRBs appears to decrease as a function 
of time (e.g. StarUng et al._2005; Gendre et al.|2007 Campana 



et al. 2007|l. For this reason, it is often preferable to choose the 
later PC data, as it is likely to be closer to the Galactic value. 
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However, the WT data often has considerably higher signal. 
Therefore, where the PC and WT mode data gave results consis- 
tent within Icr (68% confidence), both values were used. Where 
the results were discrepant at > Icr, the PC value was used. Since 
the adopted method depends solely on the limit of the population 
in the A^Hx-^v plane, the result is, however, not very sensitive to 
this. 

The cut-off, which is the relation between the X-ray column 
density and the Hi or Ay columns, was obtained from the data 
using a two-dimensional two-sample KS (2D2SKS) test. The 2D 
distribution was fit using a log-normal for the Galactic column 
density distributions (either Ay or A^hi) and the sum of the nor- 
malised Galactic column and an extragalactic column for the 
GRB afterglows, where a log-normal was also used to repro- 
duce the extragalactic column density distribution. The parame- 
ters for the Galactic column density log-normal distribution was 
obtained by fitting the Ay or A^hi distributions directly. The pa- 
rameters of the extragalactic log-normal distribution as well as 
the normalisation of the Galactic component of the GRB after- 
glow absorption distribution were left as free parameters in the 
2D2SKS fit. This normalisation parameter is the ratio between 
the X-ray absorbing column density and the Ay or A^hi- The fit 
was obtained by applying the 2D2SKS test at each step while 
stepping through each parameter across the search space, until a 
maximum in the probability was obtained. Uncertainties in the 
parameters were obtained using a Monte Carlo technique based 
on a Gaussian distribution of the errors on each datapoint. The 
same method was employed to obtain the best-fit Nui-Nux and 
Ay-Niix relations. 

3. Results 



The 
10' 



best-fit to the Av-A^Hv relation is A^h 



'^'cm"^Av 



v-jvhx reiauon is jvh = -^-^-oi 
y. The best-fit for the ratio of X-ray to 21 cm ab- 
sorbing column densities is Nux/Nhi - l-l'^o i- 

The results obtained are, not surprisingly, somewhat sensi- 
tive to finding a function to reproduce the distributions with rea- 
sonable fidelity. However, whether a normal or log-normal dis- 
tribution is used to reproduce the extragalactic column density 
distribution does not affect the results significantly. This is be- 
cause the fits are sensitive to where the limit of the 2D dataset is, 
rather than the precise shape of the distribution. As long as the 
distribution as a whole is reproduced with reasonable accuracy, 
the nature of the function is not very important. 
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In all of these fits the solar metallicity estimate of Anders 
|& Grevesse| ( |1989 | AG89) were used. It is now generally ac- 
cepted that this metallicity is ~ 45% higher than indicated by 



direct measurements of the solar spectrum ( Asplund et al. 2009 ) 
and possibly the local ISM (WAMOO). It was suggested by 
WAMOO that X-ray absorption studies should use among other 
things, updated absorption cross-sections and, in particular, up- 
dated ISM metallicity values to determine gas column densities 
in the Galaxy. The other improvements suggested by WAMOO 
have a relatively small effect on the determined gas column den- 
sities, compared to the change in metallicity they propose. Their 
proposed ISM abundances typically increase the equivalent hy- 
drogen column density for a given observation by roughly the 
inverse of the change in the metallicity. We fitted all of our 
data again using the proposed absorption model of WAMOO, i.e. 
tbabs, with abundances set at the levels suggested in that pa- 
per for the ISM. As expected, we obtain values of the Nux/Ay 
and Nhx/Nui ratios that are ~ 40% higher than using AG89 
metallicities with similar fractional uncertainties. However, pre- 
vious works have assumed metallicities similar to the AG89 val- 



Nh^ (cm-^) 

Fig. 1. The Galactic dust column (Ay) plotted against the total 
equivalent hydrogen column density measured from soft X-ray 
absorption toward Swift GRBs (Nhx)- The best-fit limit to the 
population is plotted as a solid line with Icr uncertainties as 
dashed lines and represents the Galactic relation between met- 
als and dust: Nux/Ay = 2.2;^y3 x 10^' cm^^mag"', assuming 
AG89 abundances to convert the metal column to an equivalent 
hydrogen column density. 



ues. Therefore it is the results using AG89 metallicities that we 
will use to compare with previous measurements. Indeed, as dis- 
cussed further below, it appears that the metallicity of a typical 
Galactic sightline is not consistent with the values proposed in 
WAMOO. 



4. Discussion 

The results presented here are consistent with the range of previ- 
ous results within uncertainties. Early work focussing on bright 
X-ray sources, comparing soft-X-ray metal column densities 
to exti nctio n columns by |Reina & Tarenghi| ( |1973| l, |Gorenstein| 
( [T975] l, and jRyter et ai:] ( |1975| ) pr ovided N»^/Ay = 1.9, 2.2, 



and 2.2 x 10-" cm^^mag '. Later, Predehl & Schmitt 



( 119951 



analysed 25 point sources and four supernova remnants (SNRs) 
with data from ROSAT, deriving the soft X-ray absorption from 
the spectra of the sources. That work resulted in a value of 
Nhx/Av = 1.79 + 0.03 X 10^' cm'^mag '. To extend the sam- 
ple to higher extinction lines of sight, |Vuong et al. ( 2003[ l 
analysed the X-ray absorption and 7-band extinction toward 
six nearby star-forming regions (dominated by observations of 
pOph) using pre-main sequence stars. They obtain a relation 
A^HxMy = 5.6 + 0.4 x 10^' cm^^mag"'. Converting this to 
Nux/Ay using the mean Galactic extinction curve (e.g. CardeUi 
|et al.|[T989| ) with an Ry - 3.1, yields an extremely low value, 
A^HxMv ~ 1 X 10^^ cm"^mag"'. In that paper the unusually 
low metals-to-dust ratio is ascribed to either low metallicity 
in the region of the Galaxy local to the Sun, or to a very 
flat extinction curve in the pOph cloud. The former hypothe- 
sis seems implausible since it requires that the gas-to-dust ra- 
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Fig. 2. The Galactic H i column density plotted against the to- 
tal equivalent hydrogen column density measured from soft X- 
ray absorption toward Swift GRBs (Nhx)- While the soft X- 
ray absorption is presented in units of equivalent hydrogen col- 
umn density, it is effectively a metal column density converted 
to equivalent column density assuming AG89 abundances (see 
text). The best-fit relation is super-solar and is plotted as in 
Fig. [T] above. The best-fit relation, corresponding to a ratio of 
Nhx/Nui - Ll^ii'^, would be even higher using the current 



Asplund et al. ( 2009[l solar metallicities. 



line ratios (as recently done for GRB host galaxies, |Wiersem^ 
2011') it would significantly improve the understanding of the 
extinction properties of the dust columns along the line of sight. 

The method used in this paper possesses several advantages 
over previous studies. First, the underlying X-ray source spec- 
trum is well-defined, almost invariably dominated by a power- 
law. Second, the sight-lines are passed at random through the 
Galaxy. Third, the sight-lines are to objects outside the Galaxy, 
allowing us to include the entire Galactic column in any given 
direction, and is therefore a good census of the whole Galaxy. 
However, the drawback is clearly that we are dealing with an 
upper-limit and therefore a large fraction of our more than 600 
sightlines do not contribute much statistical power to the con- 
straint on the metals-to-dust ratio. Partly for this reason, the un- 
certainties on the best-fit linear relation are fairly large, ~ 20%. 

A curious fact among previous measurements deriving the 
metals-to-extinction ratios from X-ray data is that the quoted sta- 
tistical errors are typically a few percent, while clearly the statis- 
tical quality of the fits are low, with very large outliers in terms 
of contributions to the fit statistic (e.g. SNR GO.O-hO.O in |Guver 
& Ozel 2009 shows a very large deviation from the fit). From 
this fact alone it is obvious that either the statistical or systematic 
uncertainties in these studies are substantially underestimated or 
that the intrinsic variation in the metals-to-extinction ratios is 
larger than the uncertainty. In addition, the studies quoted above 
reach mean values that di ffer by significantly m ore than a few 
percent, e.g. th e studies of Giiver & Ozel ( 2009 1 and Predehl & 
Schmitt ( 1995 1 differ by more than 4cr, supporting the fact that 
the quoted statistical uncertainties do not represent the actual 
scatter in the data. The study on the p Oph cloud is instructive in 
regard to deciding the origin of the scatter ( Vuong et al._2003) , 
providing what appears to be a good quality fit for a single re- 
gion in the Galaxy, and hints that perhaps it is the inherent scatter 
in the metals-to-extinction ratio that is responsible for the large 
deviation from one study to the next. 



tio remains constant while the metals-to-dust ratios and metal- 
licity are lower by a factor of two. It therefore seems consid- 
erably more plausible that the dust grains in a dense molecu- 
lar cloud are larger, and that therefore the extinction curve is 
simply flatter. Adopting Ry = 4- ( Vrba et ar]|1993^ re sults in 
Nh^/Av = 1.7±0.1x10^' cm"^ mag"' ; adopting Ry = 6 ( Vuong 
let al."2003^ gives 1.9 x 10"' with a similar error. Observations 



of 38 sightlines in the far-UV allowed [Rachford et al. (2009) to 
estimate the gas-to-dust ratio including the effects of molecu- 
larisation of the hydrogen. They found that for their sightlines, 
withE(B-y) in the range 0.17-1.08, the fraction of hydrogen in 
molecular form, ///,, was typically ~ 0.5 with none approaching 
fH2 ~ 1 ■ Their total gas-to-dust ratio was consistent with previ- 
ous measurements, Nn^l Ay - 2.15 + . 14 X 10^' cm"^mag"'. 
Finally, a recent paper by Giiver & Ozel ( 2009 1 analysed this is- 
sue very systematically by selecting 22 Galactic SNRs, deriving 
the metal absorption from the soft X-ray spectra and the extinc- 
tion from the Balmer decrement in most cases. They obtained 
A^HxMv = 2.21 ±0.09x 10^' cm"^ mag"'. In that work they sug- 
gest that while 143 SNRs exist in the Chandra, XMM-Newton, 
and Suzaku archives that permit good measurements of A^Hx^ 
only the 22 they use have reasonable dust column density mea- 
surements available. An obvious extension of that work would 
be to obtain extinction measurements for the remaining SNRs 
and increase the sample size from ~ 20 to over a hundred. A po- 
tential improvement of that technique could be to do more than 
use the HafH/3 Balmer decrement which gives a measure only of 
the optical reddening. If the technique could be extended to other 



4.1. The effect of metallicity 

The results presented here are compatible with previous obser- 
vations since as far as can be determined, previous X-ray mea- 
surements used metallicity values similar to AG89 to determine 
the equivalent hydrogen column density. The comparison made 
in this paper with the 21 cm Hi measurements indicates that the 
mean ISM metallicity in the Galaxy is approximately the AG89 
value, rather than the current best estimate of the solar metallic- 
ity as proposed by WAMOO. While the uncertainties presented 
in this study are relatively large, a more telling comparison is 



with Lya studies of gas-to-dust ratios in the Galaxy (Bohlin et al. 
,1978„Whittet.l98llppras & Savage|1994| e.g.), which suggest 
a value of A^hi in the range 1.6 - 1.9 X 10^'Aycm"^. The ratio 
of the typical X-ray-derived value to the Lya-derived value is 
around 1.2, immediately indicating that the mean Galactic ISM 
metallicity is somewhat larger (approximately 20%) than the so- 
lar metallicity values of AG89. Such a result is a little surprising 
given the debate in the X-ray literature on the correct dust-to-gas 
ratio. However, it has been known for many years that the ISM 
metallicity increases toward the Galactic centre and that the ISM 
metallicity becomes similar to the AG89 solar value at around 
6kpc from the Galactic centre (I Rolleston et al.|2000| . Since the 
ratios derived in X-rays are dominated by high-extinction sight- 
lines, it is perhaps after all not that surprising that the metallicity 
obtained for a typical sightline through the Galaxy is ~ 75%, or 
0.25 dex, higher than the most recent measurements of the solar 
metallicity ( |Asplund et al.|2009] l. This value is particularly use- 
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ful since it is a direct measure of the ISM itself without having to '""' ' 

use stellar photospheres and it is effectively independent of dust 22 

depletion or ionisation in the ISM. 

4.2. Proxies for the soft X-ray absorbing coiumn density 

For a typical X-ray sightline through the Galaxy, therefore, a 
more accurate estimate of the soft X-ray absorption may be ob- 
tained by using a metallicity value somewhat higher than that of 
AG89. Certainly it seems that using the current solar metallici- 
ties ( Asplund et al.|2009) l with Hi column densities is likely to 
lead to a large underestimate of the Galactic soft X-ray absorp- 
tion. 

Since the Galaxy is known to show a strong metallicity gra- 
dient fRolleston et al. 2000 1, if the dust-to-metals ratio is roughly 



constant ( Draine 2003 1, one would expect the dust-to-gas ra- 
tio to vary in a similar way, i.e. a higher dust-to-gas ratio to- 
ward the Galactic centre. Sightlines toward the Galactic cen- 
tre will have higher column densities: we therefore may expect 
a relationship in the Galaxy between the dust-to-gas ratio and 
the total column density. Such a relationship is possibly indi- 
cated by Schlegel et al. ( 1998| l, where the Galactic plane has a 
significantly higher dust-to-gas ratio than at high-latitudes with 
the Galactic centre direction being especially high. In Fig. [3] 
the relationship between the Galactic 21 cm Hi and dust col- 
umn densities is shown. Sightlines toward the Galactic centre 
{cos(l) cos{b) > 0.9, open squares on Fig.lSll do indeed seem to 
show a higher dust-to-gas ratio than sightlines with comparable 
column densities along other sightlines. The comparison of dust 
to gas directly yields a non-linear relation, with higher column 
density sightlines showing a higher dust-to-gas ratio. The linear 
relation A^hi - 2.2 x lO^'Ay shows clear residuals. However, 
A^Hi = 2.00 + 0.14 X 102'A^^^±°°^ provides a better fit to the 
data, though the fit is still not good. 



Previous studies have shown a similar effect. Schlegel et al. 
( |1998| ) note that the dust content of low density high velocity H 1 
clouds is lower than in other regions of the Galaxy. The survey 
of H I column densities and extinction in the UV using Lya ab- 
sorption and reddening of stars of Diplas & Savage] fl 994 ) show 
the same eff'ect very clearly (their Fig. 4a and to a lesser extent 
4b). 

In principle, a change in the absolute to selective extinc- 
tion ratio, Ry, could be responsible for the apparent change in 
gas-to-dust ratio derived from reddening. However, since we ob- 
serve the same effect in Fig. [3] where the extinction data are de- 
rived from the dust emission in the far-infrared (FIR), not from 
reddening, this seems unlikely. Another possibility is that dust 
temperatures, which strongly affect the emitted flux in the FIR, 
might produce such an effect - cooler lines of sight would yield 
a spuriously lower dust column if the assumed temperature was 
too high. However, [Schlegel et al.| ( [T998[ ) accounted for the dust 
temperature in their analysis, and more compellingly, we ob- 
serve the same eff'ect in extinction in the UV data of [Diplas &| 
|Savage| ( fT994l l. 

Other explanations for the correlation of gas-to-dust ratio 
with column density, such as ionisation of the hydrogen at 
low column densities or molecularisation at high column densi- 
ties, seems unlikely since the gas-to-metals ratio remains effec- 
tively constant across the range we study using H i observations 
(Fig.l2]l. It is worth noting, however, that the deviation observed 
at the high column end of Fig [3] is consistent with the factor 
of two decrease in the H i column density observed in far-UV 
molecularisation studies, where a approximately half of the hy- 
drogen exists in molecular form along sightlines with moderate 
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Fig. 3. Atomic hydrogen column density plotted as a function of 
extinction for the sightlines associated with every GRB in the 
sample from the LAB Hi survey ( Kalberla et al.|2005) and the 
dust maps of [Schlegel et al. (1998i. The anticipated linear re- 
lation is plotted as a solid line with residuals to this relation in 
terms of Ax in the lower panel. The linear relation is clearly a 
poor fit to the data, suggesting that the dust-to-gas ratio increases 
as a function of column density. The best-fit non-linear relation 
is plotted as a dashed line. The non-linearity in the relation may 
be due to metallicity variations along different sightlines. The 
metallicity gradient in the Galaxy suggests that sightlines toward 
the Galactic centre should be more metal rich and therefore show 
a higher dust-to-gas ratio. Sightlines toward the Galactic centre 
(cos(0 cos(b) > 0.9) are plotted with open squares and do indeed 
show a higher dust-to-gas ratio than sightlines with comparable 
column densities. 



or greater extinctions (Bohlin et al. 1978 ; Rachford et al. 2009J. 
Even if molecularisation or ionisation were partly responsible, 
however, it does not alter the fact that the H i column density is 
not the best proxy for the soft X-ray absorbing column density. 

This result immediately suggests that while using A^hi as a 
proxy for the soft X-ray absorption is reasonable, using the dust 
column, as suggested by Schlegel et al. ( [T998| ), is likely to yield 
better results not simply because of the higher resolution of the 
dust maps, but because the dust-to-metals ratio seems to be more 
constant than the dust-to-gas ratio. 

While this cannot be investigated further here due to the 
large uncertainty in the ratios we derive, future studies should 
try to determine the systematic Galactic radial and column den- 
sity dependences of the dust-to-metals ratio. This latter point 
may be important since it is well-known that the depletion of 
metals out of the gas phase increases as we move into the disk 
and into cooler environments (Savage & Sembach'1996), which 
would obviously be associated with higher column dens ities. It 
has been suggested above that the exploration of [Giiver & Ozel 
(2009, 1 could be extended by increasing the sample size. Such 
an extension would indeed be worthwhile and allow the Galaxy 
to be divided into various lines of sight, related to the disk or 
the bulge. A more complete approach might be to use bright ex- 
tragalactic objects known to have low host galaxy absorptions, 
which would allow a census to be taken of sightlines through 



D. Watson: The Galactic dust-to-metals ratio using GRBs 



arbitrary directions in the galaxy, including the halo. An obvi- 
ous candidate are blazars, where the dust extinction could be 
derived either as done here from Galactic dust maps, or from 
reddening of the optical power-law itself (where the caveats 
are that the data in different bands be taken simultaneously and 
that one needs to account for intrinsic curvature of the spectra). 
Indeed, for a smaller sample, the X-ray and some of the optical- 
UV data already exists. Blazars might be bright enough for UV 
spectroscopy as well, allowing H2 and Lya measurements to be 
made. Such a programme, while observationally intensive for a 
large enough sample to probe a significant fraction of the Galaxy, 
could represent one of the best probes of the Galactic ISM yet 
devised. 



4.3. Extragalactic metals-to-dust ratios 

Studies of metals-to-dust and gas-to-dust also exist outside 
the local group. Metals-to-dust ratios from foreground lensing 
galaxies at z < 1 have been obtained using multiply-imaged 
quasars ( |Dai et al.|2006[|Dai & K ochanek'2009). Results show 
metals-to-dust ratios consistent with Galactic values. Those ob- 



jects are typically relatively high mass, evolved systems. Zafar 
et al. ( 2011| l showed in an analysis of GRB afterglows that GRB 



host galaxies, which are typically young and star-forming with 
low metallicities and hard radiation envi ronments (jC hristensen 
eTaLlpOCM, Castro Ceron et al.|[2006l [Savagho et al. 2009 i 



Svensson et al.||2010[ [Chary et al.||2002[ |Fynbo etaL||20()8 
Watson et al.||2010 1, have substantially lower dust-to-gas ratios 



than the local group even after accounting for metallicity. 

5. Conclusions 

An analysis of the dust-to-metals ratio and metallicity of sight- 
lines through the Galaxy has been presented. The Galactic metal 
column densities were determined using the lower bound of 
the distribution of soft X-ray absorptions of the afterglows of 
a large sample of GRBs detected by the Swift satellite. The cor- 
responding extinction and gas column densities were found us- 
ing the dust and H i maps of Schlegel et al. ( 1998| l and Kalberla 
et al.| (12005 ) respectively. The metal to atomic hydrogen re- 
lation is well reproduced with a metallicity ~ 1.75 times the 
solar metallicity of Asplund et al. (200911. The best-fitting re- 
lation between metal and dust column densities is NhxI^v - 
2.2^Q3 X 10^' cm^^mag"' (using AG89 abundances). Previous 
observations are consistent with this result, suggesting that the 
metallicity for a typical ISM sightline is 0.25 dex higher than the 
current best value for the solar metallicity. It is therefore sug- 
gested that a better reproduction of the Galactic soft X-ray ab- 
sorption will be provided with a metallicity ~ 20% higher than 
the solar metallicity of AG89. 

However, it is also found that a linear representation does 
not reproduce the gas-to-dust relationship very well. A gas-to- 
dust relationship with A^hi - 2.00 ± 0.14 x A^^^*°'" provides 
a much better fit to the data. It is very likely that this is pre- 
dominantly a metallicity-gradient effect, and it is therefore con- 
cluded that while the gas-to-dust relation may be as given above, 
the best proxy for the Galactic soft X-ray absorption should be 
given by the dust column density with a relation of Nu^jAy - 



2.2x10^1 cm 
for an Ry = 3.1 



^ mag ' or A^i 



Hx 



I E{B-V) ^ e.'i.XlQ'^^ cmr^ mag- 
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